Singular field response and singular screening of vacancies in antiferromagnets 
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For isolated vacancies in ordered local-moment antiferromagnets we show that the magnetic-field 
linear-response limit is generically singular: The magnetic moment associated with a vacancy in 
zero field is different from that in a finite field h in the limit h — > + . The origin is a universal and 
singular screening cloud, which moreover leads to perfect screening as h — > + for magnets which 
display spin-flop bulk states in the weak-field limit. 



Defects arc ubiquitous in solids. In magnets with local- 
ized spin moments, typical classes of defects are missing 
or extra spins, arising, e.g., from substitutional disor- 
der. Very often, even small concentrations of such defects 
produce a large magnetic response at low temperatures: 
Quasi-free spins cause a Curie tail in the magnetic sus- 
ceptibility, which then is routinely subtracted from raw 
experimental data. Assuming independent defects, the 
amplitude of the Curie tail can be utilized to estimate 
the defect concentration, provided that the behavior of a 
single defect is known. 

Here we discuss the physics of isolated vacancies in 
antiferromagnets (AF) which display semiclassical long- 
range order (LRO) in the ground state [ljj. In zero mag- 
netic field, the state with a single vacancy has a finite uni- 
form magnetic moment, too, because the vacancy breaks 
the balance between the sublattices. For collincar mag- 
nets, mo is quantized to the bulk spin value, too = S 
@j 9 j while in the non-collinear case fractional values of 
too occur due to the local relief of frustration Q ■ These 
vacancy moments are expected to show up in magnetiza- 
tion measurements, and they produce a low-temperature 
Curie contribution to the uniform susceptibility in the 
two-dimensional (2d) case where bulk order is prohibited 
by the Mermin- Wagner theorem 

In this paper we show that, in an applied field h, non- 
trivial screening of the vacancy moment occurs, such that 
the linear-response limit h — y + is singular for a mag- 
net with a single vacancy, Fig. [TJ The vacancy-induced 
magnetization jumps discontinuously from its zero-field 
value too to a different value m{h — y + ) upon apply- 
ing an infinitesimal field h. Thus, measurements of the 
vacancy-induced moment m(h) in a finite field h can- 
not detect the zero-field value to even for small h 0, 
which is of obvious relevance for any experiment trying 
to quantify the defect contribution to a sample's magne- 
tization or susceptibility. Furthermore, the spin texture 
around the vacancy at finite h has a piece [9( which is 
singular as h — > + - in a sense made precise below - 
which screens the vacancy-induced moment perpendicu- 
lar to h. For magnets which feature spin-flop states (with 
all spins perpendicular to h as h —> + ) in the absence of 
the vacancy, this leads to a semiclassical version of perfect 
screening of the vacancy moment, m(h — y + ) = 0. 



In the body of paper, we present general arguments 
and microscopic calculations supporting these claims. 
Explicit results will be given in a 1/S expansion for spin- 
s' AFs on 2d lattices, with 



(1) 



but our results are valid for AFs with LRO in any di- 
mension d. In particular, the singular response occurs 
for vacancies in the square-lattice AF, where - despite 
numerous studies d, 0, [SHlH - it has been overlooked 
to date [HI]- The singular behavior will be cut off for a 
finite vacancy density, and we shall discuss the resulting 
crossover scales [1| . 

Bulk behavior in a field. To set the stage, we reca- 
pitulate the evolution of a bulk AF state, which spon- 
taneously breaks the underlying SU(2) symmetry, upon 
application of a uniform field. As h ^ breaks the sym- 
metry down to U(l), an infinitesimal field typically se- 
lects a subset of states, which we refer to as h — >• + bulk 
states: For the square lattice, these are spin-flop states 
with staggered spin directions perpendicular to h. With 
increasing field the spins rotate toward the field direction, 
until a fully polarized state is reached, see the illustration 
in Fig. Hi. 








FIG. 1: Illustration of singular screening for the square- 
lattice AF: At zero field, the vacancy induces a moment mo 
with mo = S, which is modified into rh(h) in the presence of 
a uniform field h. The size of the screening cloud, lh, diverges 
as h — > + , and m — >• in this limit. For details see text. 
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FIG. 2: Vacancy contribution to the uniform magnetiza- 
tion, m(h), for the square-lattice AF. (a) Crosses: classical 
result m(h)/S, obtained on lattices with L < 768 and extrap- 
olated to L — > oo; the circle shows the linear-response value 
Imo/S 1 ! = 1. Open dots: 1/5 correction m^°\h), here L < 50. 
Vertical dashed line: saturation field h sat = SSJ Top ar- 
rows: evolution of the bulk spin state with increasing h. Inset: 
small-ft behavior of m(h)/S, together with a fit m oc hlnh 
(see text), (b) Finite-size data m(h)/S as function of 1/L. m 
saturates for L > 1^ with lh oc SJ/h. 



For the triangular-lattice AF, an anomalously large 
ground-state degeneracy exists at the classical level for 
h ^ 0, which is lifted both by fluctuations and by addi- 
tional interactions EMU. Two cases are important: (a) 
"coplanar" and (b) "umbrella" states, see Fig. [3] In (a) 
the spins are oriented in a plane in field direction, and the 
evolution is via an intermediate up-up-down magnetiza- 
tion plateau (with a bulk magnetization of 1 /3 of the sat- 
uration value). In contrast, in (b) the spins are arranged 
in planar spin-flop states perpendicular to the field for 
h — y + , and evolve in a non-coplanar fashion continu- 
ously towards full polarization. Below, we shall employ 
an additional biquadratic exchange J2(ij) K(Si-Sj) 2 with 
relative strength k = KS 2 /J to select between the two 
cases in a classical calculation: k > (k < 0) favors 
umbrella (coplanar) states, with |fc| < 2/9 required to 
preserve the familiar 120° order at zero field |18j . 

Numerical results. We now present numerical results 
for the single- vacancy finite-field ground state of square 
and triangular AFs. We consider the vacancy contri- 
bution to the magnetization, m(/i), defined as the dif- 
ference between the total magnetizations of the system 
with and without vacancy. Figs. [2] and [3] display m(h) 
for the classical AF on the square and triangular lattices 
and demonstrate our main results: (i) m(h — > + ) does 
not reach the zero-field moment |too| m any of the cases, 
i.e., the impurity magnetization jumps upon application 
of an infinitesimal field, (ii) For the square lattice we 
find m(h — > + ) — > 0; the same happens for umbrella 
states in the triangular lattice. In contrast, m(h — > + ) 
tends to a finite value in the coplanar triangular- AF case. 
Fig. [3h also shows the next-to-leading term in a 1/5 ex- 
pansion for m(h), indicating that quantum corrections 
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FIG. 3: Vacancy-induced magnetization as in Fig. [2] but 
for the classical triangular-lattice AF with biquadratic ex- 
change of strength (a) k = —0.07, leading to coplanar states 
in a field, and (b) k = 0.1, leading to umbrella states. Here 
ftsat = 957(1 + 2/c), the shaded area in (a) indicates the mag- 
netization plateau. Insets: small-/i behavior. Note that mo 
takes a fc-dependent fractional value [J. 



do not qualitatively change these results. 

In the remainder of the paper we explain the physics 
behind these striking observations. As a first step, we 
show that the spin configurations are necessarily different 
for h — and h —> + in the presence of a vacancy. 

Vacancy: Zero-field state. The vacancy breaks the bal- 
ance between the sublattices of the host AF and locally 
distorts the bulk state. If the AF is collinear (and re- 
mains collinear upon introducing the vacancy) there is no 
distortion in the classical limit, and the vacancy-induced 
magnetization is mo = S. Quantum fluctuations arise 
from the action of St S7 terms in the Hamiltonian and 
hence conserve total spin, such that mo remains locked to 
S; however, the amplitudes of the (Si) are modified near 
the vacancy. In non-collinear AFs, the directions of the 
(Si) re-adjust in response to the vacancy because frus- 
tration is locally relieved. Then, mo takes a fractional 
value which depends on S and microscopic details Q, 
i.e., the vacancy spin gets partially screened, and both 
undercompensation and overcompensation are possible. 

In all cases, the local distortions decay algebraically 
with the distance r to the vacancy due to the presence 
of Goldstone modes. Typically this decay follows l/r d+1 
[3, El, sucn that the total (i.e. integrated) vacancy con- 
tribution to any observable is bounded, i.e., of order unity 
- this we refer to as a "regular" screening cloud. 

Singular h — > + limit. Now consider applying an in- 
finitesimal field h. To this end, rotate the h = state in 
spin space such that the spins far away from the vacancy 
match a h — >• + bulk configuration. This generates a 
state with finite magnetization mo which, however, is 
in general not oriented in field direction. This is illus- 
trated in Fig. Q] for the square lattice, where the spins 
are oriented perpendicular to the infinitesimal field - the 
same applies to mo- Clearly, such a state cannot be the 
ground state for h ^= 0, as stability demands m \\ h in the 
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presence of SU(2) symmetry. In other words, the states 
at h = and h — > + are required to differ because the 
latter is subject to the additional constraint that the uni- 
form moment has to point in field direction (while both 
states have matching spin directions far away from the 
vacancy). Fulfilling this constraint requires an additional 
distortion to screen the moment perpendicular to h. 

This analysis also reveals exceptions without zero-field 
singularity in the presence of a vacancy: This happens 
when the properly rotated zero-field state has a moment 
mo which points along the direction of the infinitesimal h. 
A concrete example is the triangular lattice with copla- 
nar states and undercompensated impurity moments 
(In the classical limit, there is overcompensation for all 
k < Q which leads to results as in Fig. [3^, but under- 
compensation can be expected for S < oo.) 

We note that mean-field theory, neglecting spatial vari- 
ations and thus distributing the screening process over 
the entire system, qualitatively reproduces the zero-field 
singularity in m(h), but fails in other respects plj . 

Finite-field screening cloud. We turn to the finite-field 
screening, in order to understand the difference between 
the cases in Figs. [2^ and[3k, where m(h — > + ) reaches 
zero or a finite value, respectively. The vacancy-induced 
modification of the directions of (Si) can be parameter- 
ized by the spherical angles <50j, 6<f>i (one angle is suf- 
ficient for the coplanar cases in Figs. [2jt andEK). Con- 
sidering that the distortion serves to screen the vacancy 
moment's component perpendicular to h, the screening 
can be quantified using 

m±(r) = J2 (Si) ± (2) 

\ri I <r 

where (Si)± is the moment perpendicular to h, and the 
sum runs over all sites with a distance to the vacancy 
closer or equal r. Clearly, m±(r — > oo) is the total mag- 
netization perpendicular to h and has to vanish. 

We start with the square lattice: The angles <50, 
Fig. [4ji, decay exponentially on a length scale lh oc SJ/h. 
This is expected: The distortion is mediated by a bulk 
mode which is a gaplcss Goldstone mode in zero field, 
but acquires a gap oc h for h > 0. A remarkable fea- 
ture is that 5Qi at fixed r, has a non-monotonic field 
dependence and appears to reach zero as h — > (visible 
at small r). This implies that the distortion, and with 
it the screening, happens at progressively larger r with 
decreasing h. This is confirmed by the plot of m±(r), 
Fig-Hti showing that the cloud size increases as h — > + . 

To understand these results analytically, we resort to 
a continuum description: The small-field distortions con- 
sist of smooth variations of the order-parameter field <p(r) 
(or, equivalently, of the angles on all sublattices). For 
a single vacancy, those distortions have been analyzed 
in Ref. 0, and found to decay as h d ~ Y / '(p s c d ~ 2 ) fd(hr / c) 
where p s and c oc SJ are the spin stiffness and spin- 
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FIG. 4: Field-induced screening cloud for (a,c) the square 
lattice and (b,d) the triangular lattice with k — —0.07, for 
different values of h. (a,b) Vacancy-induced rotation angles 
SQ(f) as function of the distance r to the vacancy. In (b) the 
zero-field angles are shown as well. Insets: scaling of 50 /h 
vs. rh/(SJ) oc r/lh, with the short-range piece projected out 
in (b). (c,d) Transverse magnetization m±(r) as function of 
the integration radius r in Eq. ([2]). Insets: fn± vs. rescaled 
distance rh/(SJ). 



wave velocity, respectively, and fd is a universal func- 
tion which depends on the dimension d only. Specifi- 
cally, fz(x) = Kq(x) (the modified Besscl function), and 
f 3 (x) = exp(-x)/x. 

To analyze the screening of m±, we first note that 
distortion-induced (local) contributions to m±_ arise - 
in the presence of a finite bulk magnetization density 
I oc h - to linear order in 50, 5m± oc J d d r\l\8Q, which 
have to compensate the bare vacancy-induced trans- 
verse moment (which itself equals S in the square-lattice 
case). Using the above long-distance form of the dis- 
tortion, we have Sm± oc (c 2 / p s )(h/c) d J d d rfd(hr/c) oc 
(c 2 /p s ) J d d xfd{x). As the last integral is finite, we see 
that all h dependence has dropped out, indicating stable 
screening of m± for any h. 

The insets in Figs.0k,,c confirm the anticipated scaling 
of our data with 1^. Plotting SO/h and m± each as func- 
tion of rescaled distance rh/(SJ) reveals data collapse, 
i.e., the screening cloud is universal for small fields. The 
scaling of SO /h also implies that the texture does not 
approach a well-defined limiting form as h — > + , but 
instead has a diverging size combined with a vanishing 
amplitude - this we refer to as "singular" screening cloud. 

A similar calculation can be employed to obtain 
the cloud's contribution to the field-parallel magnetiza- 
tion mil = m(h). This arises from two higher-order 

terms, J d d r\l\SQ 2 and J d d r\(p\V 2 0, while lower or- 
ders vanish because the transverse magnetization den- 
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sity is zero in the bulk reference state. The first term 
scales as h d_1 J d d xf%{x) while the second one scales as 
h J d d xf'[(x). Importantly, the first integral is regular 
whereas the second one has a short-distance logarithmic 
divergence in any d - this is cut off by noting that the 
lower limit of the x integral is not zero, but oc h. As a 
result, the cloud's contribution to m{h) follows (h\nh) 
for small h. Indeed, our numerical results in Fig. [2t are 
perfectly fit by this form - note that the bare vacancy- 
induced moment along h vanishes oc h and is subleading. 

Now we turn to the triangular lattice. An analysis of 
the angles 8<d (r), Fig. HJd, shows that they arise from a 
superposition of two textures, with spatial s-wave and 
/-wave symmetries, respectively. The latter is well lo- 
calized near the vacancy and evolves smoothly into the 
zero- field texture [H, while the former is a field-induced 
distortion very similar to that in the square lattice. As a 
result, we find two-stage screening at small h: At short 
distances, the moment of the missing spin is screened 
to mo in a regular and weakly ficld-dcpcndcnt fashion, 
while the field-perpendicular component of this too is 
then compensated at larger distances via a singular cloud 
of size lh- This singular piece can be isolated by project- 
ing SQ(r) onto its s-wave component, Fig. Hp, and it ad- 
mits the same continuum description as above. Because 
of the nearby magnetization plateau, universal scaling is 
reached at much smaller h as compared to the square 
lattice. 

Finally, this insight allows to deduce the value m(h — !> 
+ ): This is simply the field-parallel part of mo in the 
(properly rotated) zero-field state, because the contribu- 
tion of the singular screening cloud vanishes as (h\nh). 
For the triangular lattice with overcompensated impu- 
rity, this too forms a 60° angle with h [3| such that 
m(h — > + ) = |?7io|/2, matching our numerical result 
in Fig. [3^, while for all systems with spin-flop states 
m(h — > + ) =0, i.e., the vacancy moment is perfectly 
screened. 

Quantum effects. Although the concrete calculations 
so far were for classical spins, the results qualitatively 
apply to any system with LRO: As happens in general 
for non-collinear magnets, quantum effects will modify 
directions and amplitudes of {Si). However, the exis- 
tence of both the zero-field singularity and the perfect 
screening are unaffected, as they derive from symmetry 
and stability arguments. 

To illustrate this, we have calculated quantum correc- 
tions to m{h) in a 1/5 expansion using spin-wave theory 



on finite lattices [20|, [22|, [23|. A sample result for the 
square-lattice S° term of m{h) is shown in Fig. [2] [l9| - 
its ft. — >• behavior is consistent with {h In h) . 

Limits and crossover scales. The zero-field singularity 
only occurs in the thermodynamic limit to evolves 
smoothly for a single vacancy in a finite-size system, and 
the limits N — > oo and h — > do not commute 211. 



singularity is replaced by a crossover governed by the two 
length scales lh and l lmp = n i J L ^ d , the mean vacancy dis- 
tance. For lh <C 'imp, the vacancy moments respond inde- 
pendently to the field, whereas for lh 3> 'i mp the screening 
clouds overlap and hence their moments tend to average 
out (for an equal distribution over all sublattices). Thus, 
the magnetization per vacancy will follow m(h) as calcu- 
lated here for elevated fields, but will (smoothly) vanish 
below a crossover field given by h/ J ~ ^j^p- 

These considerations also apply at finite temperatures 
inside the ordered phase of 3d Hcisenberg magnets. In 
2d, order occurs only at T = 0, but the bulk correla- 
tion length £ becomes exponentially large as T — >• 0. 
Then, the physics is governed by the interplay of the 
three length scales lh, 'imp, and £, and we quickly discuss 
some interesting limits. For £ 3> lh, 'imp one recovers the 
T = physics discussed above. Single-impurity physics 
obtains in the dilute limit, 'i mp 3> £, lh, where at elevated 
fields, £ 3> lh, the magnetization per vacancy again fol- 
lows our m{h). In contrast, in the low-field limit, £ <C lh, 
linear response is restored such that the vacancy suscep- 
tibility takes the Curie form x = rn%/(3kT) i.e., 
the singularity is replaced by a crossover on the scale 
h/J ~ cxp(— T/J). Finally, the limit lh £, 'i mp is gov- 
erned by zero-field physics. 

Conclusions. For a single vacancy in an ordered AF we 
have found that the magnetic behavior is generically sin- 
gular in the weak-field limit: The vacancy contribution to 
the magnetization jumps discontinuously upon applying 
an infinitesimal field. The singularity can be traced back 
to non-trivial field-induced screening due to a universal 
and singular screening cloud in the h — > + limit. 

Our predictions can be experimentally verified in any 
AF doped with a small, controlled amount of vacancies, 
provided that magnetic anisotropies are small Q . In ad- 
dition, numerical studies beyond the 1/5* expansion (e.g. 
quantum Monte Carlo) are called for. 

We thank J. Kunes, R. Moessner, and O. I. Motrunich 
for discussions, and L. Fritz for a collaboration at an early 
stage of this work. This research was supported by the 
DFG through FOR 960 and GRK 1621. 



Assuming now a finite vacancy concentration n; m p , the 
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I. VACANCY IN A FINITE-SIZE SYSTEM 



II. MEAN-FIELD APPROXIMATION FOR A 
VACANCY-DOPED MAGNET 



While the analysis in the main text was concerned with 
the vacancy magnetism in an infinitely large host, we 
discuss here the behavior for a finite host of size N = L d 
where L is the linear dimension, i.e., where the vacancy- 
doped system has (TV — 1) spins. 

The zero-field singularity advertised in the paper is 
due to competing boundary (i.e. impurity) and bulk 
magnetization processes, described by terms in the sys- 
tem's energy of the form — hm and — xh 2 N , respectively, 
where x is the bulk susceptibility. Importantly, the lim- 
its N — > oo and h — > do not commute: If h is taken to 
zero in a finite-size system, then the impurity contribu- 
tion to the field energy will always dominate over the bulk 
contribution for sufficiently small h, such that the zero- 
field vacancy moment mo simply becomes field-aligned 
for h <C SJ/N . In this limit, the bulk spins will not 
have taken their "spin-flop" configuration. Upon increas- 
ing h, the bulk term gradually starts to compete with 
the impurity term, leading to a smooth evolution of m in 
a finite-size system, with a crossover between small-field 
and large-field behavior at h ~ SJ/N. 

In contrast, in the thermodynamic limit where N — > oo 
is taken first, the bulk contribution to the field energy 
dominates for any finite h, such that the bulk spins are 
in their "spin-flop" configuration which forces m(h) to be 
different from mo as discussed in the main text. In other 
words, the above finite-size crossover in m(h) becomes 
the advertised jump occurring at h = + in the thermo- 
dynamic limit. 

These considerations are well borne out by numeri- 
cal simulations on small systems in the classical limit, 
Fig. ISU where m(h) follows an initial decrease from mo 
for small h before it sharply crosses over into its N — > oo 
behavior for larger h. For the square lattice, the crossover 
field is given by h x = 8SJ/N. 

We note that this finite-size crossover is irrelevant for 
experiments in macroscopic solids, because there any ap- 
plied field is larger than J/N. Instead, one needs to con- 
sider the crossovers induced by a finite concentration of 
vacancies, which are discussed toward the end of the main 
text. The finite-size crossover may, however, be probed 
in finite systems realized, e.g., with ultracold gases in 
optical lattices. 



For completeness, we describe here the results of a 
mean-field approach to the (classical) vacancy-doped an- 
tiferromagnet. The simplest approximation assumes all 
spins on a given sublattice to be parallel, such that each 
sublattice can be represented by a single "superspin". For 
the square lattice, this amounts to a mean-field energy of 
the form 

Emf = ^fs A ■ S B - h(S z A + S B ) (SI) 

where z = 4 is the number of nearest neighbors and Sa, b 
represent the magnetization vectors of the sublattices A 
and B, with \S A \ = SN/2 and \S B \ = S{N/2 - 1) for the 
single-vacancy case. Minimizing -Emf w.r.t. the angles 
parameterizing the directions of Sa. b yields a mean-field 
approximation, e.g., for the field-dependent magnetiza- 
tion of the vacancy-doped antiferromagnet from which 
one can obtain the vacancy-induced magnetization m(h) 
as defined before. 

Formally, Eq. (|Slj) is a mean-field theory for a 
ferrimagnetii any spatial variations which are expected 
due to the presence of the impurity are discarded. Such 
a mean-field theory has been briefly considered for the 
vacancy problem in Ref. [U- 

Mean-field results for m(h) for finite square lattices are 
included in Fig. ISlb . As compared to the numerically 
exact results, various features are apparent: 

• The mean-field results are extremely accurate for 
small fields h < 8SJ/N. This can be rationalized 
by considering that, in the full solution, the field- 
induced length scale of the spin texture, given by 
lh oc SJ/h, then becomes lh > N in this regime. 
Hence, 3> L, which implies that there is essen- 
tially no spatial variation of the spin pattern across 
the entire system, and consequently a theory ne- 
glecting spatial variations is accurate. 

• The mean-field theory incorrectly predicts m(h) = 
for all intermediate fields 8/N < h/(SJ) < 8 - 
8/N. 

• Taking the limit N — > oo, the mean-field theory 
reproduces the correct jump in m(h) at zero field, 
but incorrectly predicts a jump at h = h~ at as well. 

For the triangular lattice, the mean-field approach can 
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FIG. SI: Vacancy contribution to the uniform magnetization, m(h) (red), calculated for a classical finite system of size N — L 2 , 
for (a) the square-lattice AF, (b) the triangular-lattice AF with k — —0.07, and (c) the triangular-lattice AF with k — 0.1. In 
all cases, the left panels show the full field range, while the smaller right panels magnify the behavior near h — and h = h sa ,t- 
The corresponding infinite-system results are in Figs. 2 and 3 of the main paper; the shaded region in panel (b) indicates the 
magnetization plateau. Also shown are the results for m(h) from the mean-field approach (green) as given by Eqs. (|S1IS2[) for 
the same parameters and system sizes, for details see text. 



be generalized to three sublattices, where now 

3zJ -> -> - 

Emf = -^rr(SA-SB+SA-Sc+SB-Sc) — h(S%+S% + Sc) 

(S2) 

with z = 6 and \S AtB \ = SN/3 and \S C \ = S(N/3-l) for 
a single vacancy. An additional biquadratic interaction 
(see main text) is straightforward to include. 

Figs.lSlb.c show the mean-field results for m(h) in the 
triangular lattice with k = —0.07 and k = 0.1 as in the 
main paper, vis-a-vis the exact results. While the mean- 
field theory again predicts a jump of m{h) at h = + in 
the N —> oo limit, its deficiencies are even more apparent 



than in the square-lattice case: The mean-field theory 
fails to properly capture the non-singular piece of the 
screening process, such that both m and the jump height 
are incorrect. Further, large deviations occur near the 
magnetization plateau and the saturation field. 

In summary, the mean-field approach qualitatively re- 
produces the small-field finite-size crossover in m(h) and 
its zero-field jump in the thermodynamic limit, essen- 
tially because it is bound to follow the general arguments 
laid out in the main text. However, the neglect of spa- 
tial fluctuations renders the mean-field theory unreliable 
for various important aspects of the impurity problem at 
hand. 
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